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Abstract 

Background: Serotonin plays a pivotal role in regulating and modulating physiological and behavioral processes in both 
vertebrates and invertebrates. In the honeybee (Apis mellifera), serotonin has been implicated in division of labor, visual 
processing, and learning processes. Here, we present the cloning, heterologous expression, and detailed functional and 
pharmacological characterization of two honeybee 5-HT 2 receptors. 

Methods: Honeybee 5-HT 2 receptor cDNAs were amplified from brain cDNA. Recombinant cell lines were established 
constitutively expressing receptor variants. Pharmacological properties of the receptors were investigated by Ca 2+ imaging 
experiments. Quantitative PCR was applied to explore the expression patterns of receptor mRNAs. 

Results: The honeybee 5-HT 2 receptor class consists of two subtypes, Am5-HT 2ol and Am5-HT 2 p. Each receptor gene also 
gives rise to alternatively spliced mRNAs that possibly code for truncated receptors. Only activation of the full-length 
receptors with serotonin caused an increase in the intracellular Ca 2+ concentration. The effect was mimicked by the agonists 
5-methoxytryptamine and 8-OH-DPAT at low micromolar concentrations. Receptor activities were blocked by established 5- 
HT receptor antagonists such as clozapine, methiothepin, or mianserin. High transcript numbers were detected in exocrine 
glands suggesting that 5-HT 2 receptors participate in secretory processes in the honeybee. 

Conclusions: This study marks the first molecular and pharmacological characterization of two 5-HT 2 receptor subtypes in 
the same insect species. The results presented should facilitate further attempts to unravel central and peripheral effects of 
serotonin mediated by these receptors. 
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Introduction 

The biogenic amine serotonin (5-hydroxytryptamine, 5-HT) has 
important physiological functions including the regulation of 
energy balance and food intake, gastrointestinal and endocrine 
function, and cardiovascular and pulmonary physiology (for a 
recent review, see: [1]). In the central nervous system (CNS), 
serotonin contributes to the regulation of appetite, mood, sleep, 
and cognitive functions, including learning and memory. In 
humans, impaired serotonin signaling has been implicated in 
diseases such as migraine, schizophrenia, and depression (for 
reviews, see: [2,3]). At present, 14 distinct genes encoding 5-HT 
receptors have been identified in total in mammals. These proteins 
have been assigned to seven receptor classes (5-HT! to 5-HT 7 ; 
[4,5]). Except for 5-HT :i , which is a ligand-gated ion channel, all 
mammalian 5-HT receptors belong to the family of G protein- 
coupled receptors (GPCRs). The 5-HT! and 5-HT 5 receptors 
couple preferentially to G ;/o proteins and inhibit cAMP synthesis. 



The 5-HT 2 receptors activate G q /n proteins, which mediate the 
hydrolysis of inositol phosphates and cause a subsequent increase 
in cytosolic Ca 2+ ([Ca 24 ];). The 5-HT 4 , 5-HT 6 , and 5-HT 7 
receptors couple to G s proteins and promote cAMP formation. 
Orthologous G protein-coupled 5-HT receptors with conserved 
signaling pathways have been discovered in protostomes including 
nematodes [6], crustaceans [7,8], and insects [9-11]. To date, four 
5-HT receptor subtypes have been characterized in the fruit fly 
(Drosophila melanogaster). Phylogenetically, they cluster within the 
mammalian 5-HT, A [12], 5-HT 2 [13], and 5-HT 7 [14] subfam- 
ilies. 

In addition to the genetic model organism D. melanogaster, the 
honeybee (Apis mellifera) has gready contributed to improving the 
general understanding of serotonergic control and modulation of 
various types of behavior [11,15-19]. The distribution of 
serotonergic neurons has been mapped precisely in the CNS of 
this social insect [20,21], and the serotonin content in the brain 
has been determined at various developmental stages [22,23]. 
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Serotonin levels are not constant during a bee's lifespan, with older 
foraging bees having higher serotonin contents than younger bees 
working in the hive [23]. In single-cohort colonies, which consist of 
same-aged bees, significant differences in serotonin levels have 
been found in the antennal lobes of forager compared with nurse 
bees [24], Together, these results suggest that differences in 
serotonin levels are related to specific task(s) that bees perform, 
rather than to their age. 

Several aminergic receptors have been characterized from 
honeybees in recent years [25—31]. However, our present 
knowledge about the molecular and pharmacological properties 
of 5-HT receptors remains limited. So far, only a 5-HT 1A [32] and 
a 5-HT 7 [33] receptor have been uncovered. The aim of this study 
was the molecular identification and detailed pharmacological 
characterization of honeybee 5-HT 2 receptors. 

So far, not much is known about the functions of 5-HT 2 
receptors in insects. In D. melanogaster, the Dm5-HT 2o( receptor 
(CGI 056) regulates cellular movements during germband exten- 
sion and cuticular formation during early embryogenesis [34,35]. 
Furthermore, 5-HT 2 receptors have been implicated in the 
modulation of the daily activity pattern, anticipatory behavior, 
and aggression (for a review, see: [19]). Recently, the Cv5-HT 2o( 
receptor of the blowfly, Calliphora vicina, has been shown to be 
expressed not only in the brain but also in the salivary gland [36]. 
Here, its activation by serotonin leads, via an increase in [Ca 2+ ] ; , to 
an elevation of the Cl~ permeability of both the basolateral and 
apical membrane and thus facilitates Cl~ movement from the 
haemolymph into the lumen of the gland [36] . 

Inspection of the completely sequenced honeybee genome [37] 
has revealed the existence of two candidate genes encoding 5-HT 2 
receptors: Am5-HT 2a and Am5-HT 2 p [10,1 1]. After heterologous 
expression, both receptors cause increases in ([Ca 2+ ] ; ) upon 
stimulation with nanomolar concentrations of serotonin. These 
responses are efficiently blocked by 5-HT receptor antagonists, but 
with subtype-specific patterns of efficacy and potency. Because of 
their preferential expression in glandular tissues, both 5-HT 2 
receptor subtypes are likely candidates for the control or 
modulation of important secretory processes in the honeybee. 

Materials and Methods 

Cloning of Am5-ht2 cDNAs 

Single Apis melltfera drone brains were used to prepare poly(A) + 
RNA with the Micro-FastTrack™ 2.0 Kit (Invitrogen, Karlsruhe, 
Germany). Drones possess a haploid genome and, therefore, single 
nucleotide substitutions in cDNA clones cannot be due to allelic 
polymorphisms. Synthesis of cDNA employed the AccuScript™ 
High Fidelity 1 st Strand cDNA Synthesis Kit (Stratagene, 
Amsterdam, Netherlands). Specific primers (Table SI) allowed 
the entire coding region of the receptors to be amplified. The 
polymerase chain reaction (PCR) was carried out for 2.5 min at 
94°C (1 cycle), followed by 35 cycles of 40 s at 94°C, 40 s at 54°C 
{Am5-ht2a) or 64°C {Am5-ht2P), 150 s at 72°C, and a final 
extension of 10 min at 72°C. PCR products were cloned into 
pGEM-T vector (Promega, Mannheim, Germany) and subse- 
quently sequenced (AGOWA, Berlin, Germany). The nucleotide 
sequences of Am5-ht2a and Am5-ht2fi have been submitted to the 
European Bioinformatics Institute (EBI) database (accession nos. 
FR727107 and FR727108, respectively). 

Multiple sequence alignment and phylogenetic analysis 

Amino-acid sequences used for phylogenetic analysis were 
identified by protein-protein Basic Local Alignment Search Tool 
(BLAST) searches of the National Center for Biotechnology 



Information (NCBI) database with the deduced amino acid 
sequence of Am5-ht2a (Am5-HT 2ol ) as "bait". Values for identity 
(ID) and similarity (S) were calculated by using the BLOSUM62 
substitution matrix in BioEdit 7.0.5. MEGA 4 [38] was used to 
calculate the genetic distances between the core sequences and to 
construct maximum parsimony trees with 2000-fold bootstrap re- 
sampling. The D. melanogaster rhodopsin 1 (ninaE) and FMRFamide 
receptor sequences were used as out-groups. 

Quantitative real-time PCR 

Samples from various tissues of individual honeybee workers 
and drones were collected, immediately frozen in liquid nitrogen, 
and stored at -80°C until use. Total RNA was extracted by using 
the RNeasy Mini Kit (Qiagen, Hilden, Germany) and served as 
the template for cDNA synthesis. From each sample, two 
independent cDNA syntheses were performed by using Super- 
ScriptHI (Invitrogen) according to the manufacturer's instructions. 
Quantitative real-time PCR (qPCR) was carried out on a Rotor 
Gene Q (Qiagen) by using TaqMan technology with various 
fluorescent dyes to allow duplex measurements of receptor and 
reference gene expression. The fluorescent dyes used as 5'mod- 
ifications were 6-FAM-phosphoramidite, Cy5, and Cy5.5. Black- 
Berry quencher was attached to the 3'-end of TaqMan probes. 
The sequences of the primers and TaqMan probes are presented 
in Fig. S 1 . The PCR was performed with an initial step at 60°C for 
1 min and a denaturation step at 95°C for 5 min, followed by 45— 
55 cycles at 95°C for 20 s and at 60°C for 60 s. Tissue samples of 
at least five individual bees were examined in triplicate. Mean 
copy numbers were calculated by using Rotor Gene Q software 
(Qiagen). Samples from each tissue were tested in at least two 
independent qPCR runs. Receptor transcript levels were normal- 
ized to elongation factor la (Amef-la) transcript levels (= 100%) 
by using the standard curve method. The standards covered copy 
numbers from 10 ? - 10 7 . To detect statistically significant 
differences in gene expression, one-way ANOVA followed by 
Bonferroni's multiple comparison test was applied on log(x + 1) 
transformed data as implemented in SPSS 19 (Rel. 19.0.0, IBM 
SPSS Statistics). 

Construction of expression vectors 

Expression-ready constructs of Am5-ht2oc and Am5-ht2ocAIII 
were generated in pcDNA3.1 vector (Invitrogen). PCR was 
performed with specific primers (Table SI). PCR products were 
digested with JVM and EcoRl and subcloned. The resulting 
constructs were named pcAm.5-ht2tt.-YiA and pcAm5-ht2ocA III-HA. 

In order to generate an expression vector for Amb-ht2f>, a PCR 
with specific primers (Fig. SI) was performed. The PCR product 
was digested with Hindlll and Agel and subcloned into pcDNA6/ 
myc-His A vector (Invitrogen) yielding pcAm5-ht2 ft-Yiis. All insert 
fragments were checked by DNA sequencing. 

Functional expression of honeybee 5-HT 2 receptors 

Approximately 8 u.g of pcAm5-M2a.-Y{A, pcAm5-ht2ocA III-HA, 
and pcAm5-ht2 ft-Yiis were transfected into exponentially growing 
HEK 293 cells (~4xl0' cells per 5-cm Petri dish) by a modified 
calcium phosphate method [39]. Stably transfected cells were 
selected in the presence of the antibiotic G418 (0.8 mg/ml; 
pcDNA3.1 vector) or blasticidin (0.1 mg/ml; pcDNA6 vector). 
Isolated foci were propagated and analyzed for the expression of 
the receptor proteins either by immunocytochemistry and Western 
blotting (see Fig. S2) or by functional Ca 2+ imaging upon receptor 
activation. 
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Functional analysis of honeybee 5-HT 2 receptors 

The ability of the honeybee receptors to activate G q proteins 
was assessed by monitoring changes in [Ca 2+ ]i with the Ca 2+ - 
sensitive fluorescent dye Fluo-4 (Invitrogen). Non-transfected 
HEK 293 cells and cells expressing Am5-HT 2 o(, Am5-HT 2a Aiii, 
or Am5-HT 2 p were grown in 96-well plates to a density of 
~3xl0 4 per well. Cells were loaded at room temperature with 
Fluo-4 as described earlier [40] . Each 96-well plate was transferred 
into a fluorescence reader (FLUOstar Galaxy/Optima; BMG 
Labtech, Offenburg, Germany) to monitor Fluo-4 fluorescence. 
The excitation wavelength was 485 nm, and fluorescence emission 
was detected at 520 nm. Various concentrations of biogenic 
amines or receptor ligands were added, once Fluo-4 fluorescence 
had reached a stable value in each well. The changes in Fluo-4 
fluorescence were recorded automatically. Dose-response curves 
for putative agonists/antagonists were established in at least two 
independent experiments. 

Results 

Cloning and sequence analysis of 5-HT 2 receptors from A. 
mellifera 

Two candidate genes encoding putative serotonin 5-HT 2 
receptors were identified in the completely sequenced honeybee 
genome (Aml6 and Aml7; [10]). Here, we used this sequence 
information and applied a PGR-based strategy to amplify the full- 
length cDNAs. The Am5-ht2x cDNA (Am 16 in [10]) contains an 
open reading frame (ORF) of 1,962 bp and encodes a protein of 
653 amino-acid residues (Am5-HT 2a ) with a calculated molecular 
weight of 71.5 kDa. The hydropathy profile and topology 
predictor Phobius [41] suggest seven trans-membrane (TM) 
domains, which are characteristic of GPCRs. The TM segments 
are flanked by an extracellular N-terminus of 140 and an 
intracellular C-terminus of 96 residues. An alignment of the 
cDNA with the genomic database (release Amel4.5) revealed that 
the Am.5-h(2.a gene contains seven exons and six introns (Fig. 1), 
and that it is located on chromosome LG9. Am5-ht2a is the 
ortholog of the D. melanogaster 5-ht2a gene (CGI 056) with which it 
has three introns in common (Fig. S3A+C). A remarkable 
difference exists between the cloned cDNA and the annotated 
Am5-ht2a gene sequence (accession no.: XM394798). Exon V of 
the cDNA consists of 270 bp, whereas the annotated exon V 
contains only 220 bp (Fig. 1, V annot ). This discrepancy originates 
from the usage of a different 5' splice site during the in silico 
assembly of the gene, such that the annotated exon V starts at 
position 1,202 of the Am5-h&a ORF. The 3' splice sites of both the 
cloned and the annotated sequences are completely conserved 
(Fig. 1). Notably, in all PCR experiments performed, we never 
obtained a fragment as predicted by the annotation. 

However, we isolated a splice variant of Am5-ht2oc mRNA (Am5- 
hflvAIII). In Am5-ht2xAIII, exon III (225 bp) is absent, which 
results in a shortened ORF of 1,737 bp (Fig. 1). This mRNA 
variant can potentially be translated into a protein of 578 amino 
acid residues. Compared with the full-length receptor, this protein 
will lack TM3, TM4, and the entire CPL2 (see Fig. 1). 

In addition to Am5-ht2oc, three partial mRNAs that might 
encode proteins have been annotated in silico. However, with 256 
bp (XM_001 119970), 522 bp (XM_001 122856), and 1,053 bp 
(XM_624894), these sequences are too short to code for complete 
GPCRs. Nevertheless, the deduced amino-acid sequences showed 
striking sequence similarity to a 5-HT receptor of the spiny lobster, 
Panulirus interruptus (Pan5-HT 2 p; [7]; accession no.: AY550910). We 
used the crustacean sequence to design primers annealing in close 
proximity to the start and the stop codons of the hypothetical 



honeybee ortholog. A single cDNA, which contained an ORF of 
2,202 bp, was amplified on honeybee brain cDNA {Am5-ht2p). The 
ORF is distributed on three genomic contigs (Fig. 1). Contig 
NW_00 1252966 has been allocated to chromosome LG10, 
whereas the other two, i.e., NWJ301262451 and 
NW_00 1262048, have not been assigned to any chromosome, as 
yet. The coding sequence of the gene is interrupted by five introns 
(Fig. 1) with two of them being in common with the orthologous 
gene in D. melanogaster (CG42796, Fig. S3B+C). The ORF of the 
Am5-ht2fi gene encodes a protein (Am5-HT 2 p) of 733 amino-acid 
residues with a calculated molecular weight of 80.7 kDa. The 
Am5-HT 2 p receptor contains an extremely long CPL3 of 399 
residues and a relatively short C-terminus of 25 residues. 

Similar to Am5-ht2cr., we identified a splice variant of the Am5- 
ht2p mRNA. In this variant (Am5-ht2pAII), exon II is completely 
missing. This results in a frameshift of the ORF (ORF1) with an 
early stop codon at position 1,075 - 1,077 (Fig. 1). ORF1 is 
predicted to code for a protein of 358 amino-acid residues and a 
calculated molecular weight of 39.5 kDa. The sequence from Met! 
to Gly 132 including TM1, TM2, and the extracellular loop (ECL) 1 
is identical to the full-length Am5-HT 2 p protein. Because of the 
frameshift, the consecutive sequence (Gly 133 to Pro3 57 ) differs 
completely from the full-length receptor protein. A second ORF 
(ORF2, 1,605 bp) starts within exon III and lasts to the stop codon 
of the full-length mRNA. This reading frame is predicted to 
encode a protein of 491 residues with a calculated molecular 
weight of 54.1 kDa. The amino-acid sequence is completely 
identical to that of full-length Am-5-HT 2p from Met 24:i (ECL2 
prior to TM5) to Arg 733 but obviously misses a large part of the N- 
terminal region encompassing TM domains 1-4 of an intact 
GPCR (Fig. 1). 

Sequence motifs, which are essential for the three-dimensional 
structure, ligand binding, and signal transduction of the receptor, 
are well conserved in both full-length Am5-HT 2o( and Am5-HT 2 p 
proteins (Fig. 2). An exception is the Asp residue within the D-R-Y 
motif required for G protein coupling [42] ; this is substituted by 
Gly in Am5-HT 2o( (G 237 -R 231i -Y 2 39, Fig. 2). Five and four 
consensus motifs for potential N-glycosylation (N-x-[S/T]) are 
located in the extracellular N-terminus of Am5-HT 2o , and Am5- 
HT 2 p, respectively (Fig. 2). Cysteine residues in the C-terminus 
(Cys 5 7o, Cys 372 , and CyS5 94 in Am5-HT 2ot , Cys 722 in Am5-HT 2 p, 
Fig. 2) are possible sites for post-translational palmitoylation. 
Whereas no consensus site for phosphorylation by protein kinase A 
(PKA, R-[R/K]-x-[S/T]) or protein kinase C (PKC, [S/T]-x-[R/ 
K]) is found in Am5-HT 2o( , five PKA consensus sites and four 
PKC consensus sites are present within CPL3 of Am5-HT 2 p (Fig. 
2). 

A comparison of Am5-HT 2o , and Am5-HT 2 p amino-acid 
sequences with NCBI databases identified several orthologous 
protostomian and deuterostomian 5-HT 2 receptors. For Am5- 
HT 2a , the highest amino acid identity/similarity (ID/S) existed to 
the 5-HT 2 receptor ofD. melanogaster (Dm5-HT 2o( ; [13]; ID 28.4%, 
S 34.8%). For Am5-HT 2 p, homology was more pronounced to 5- 
HT 2 receptors from the crustaceans Panulirus interruptus (Pan5- 
HT 2P ; [7]; ID 31.2%, S 42.8%) and Procambarus clarkii (Pro5-HT 2 ; 
[8]; ID 33.2%, S 40.4%). In phylogenetic tree analyses, both Am5- 
HT 2a and Am5-HT 2 p always grouped with protostomian and 
deuterostomian 5-HT 2 receptors (Fig. 3). 

Tissue-specific expression patterns of the Am5-HT 2 genes 

The tissue-specific distribution of Am5-ht2oc and Am5-ht2fi 
mRNA was determined by qPCR in tissue samples from pollen 
foragers and drones (Fig. 4, Fig. S4). Application of splice-variant- 
specific primers (Fig. SI) allowed us to unravel the expression 
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Figure 1. Genomic organization of Am5-ht2a and Am5-ht2fi genes. The Apis mellifera genomic database (Amel_4.5) was screened with the 
Am5-ht2z and Am5-ht2fi cDNA sequences. The figure shows the splicing pattern of the Am5-ht2z (upper part) and Am5-ht2(i (lower part) genes. 
Nucleotide sequences of exons are given in upper-case letters and nucleotide sequences of introns are given in lower-case letters. Full-length and 
truncated receptor variants are displayed schematically. Boxes represent exons, and bent lines represent introns. The position of the seven trans- 
membrane domains is highlighted by black bars. Exons containing a frameshift are indicated in gray. The seven exons of the Am5-ht2x gene are 
completely covered in the genomic contig NW_001 253562. The six exons of the Am5-ht2(i gene are dispersed in three genomic contigs (A: 
NW_001 262451, B: NW_001 252966, und C: NW_001 262048). 
doi:1 0.1 371 /journal.pone.0082407.g001 



pattern of the different splice variants for each receptor gene. PCR 
fragments originating from full-length and splice variants of Am5- 
ht2a and Am5-ht2f} were found in the CNS, the hypopharyngeal 
glands, the salivary glands, and the Malpighian tubules of pollen 
foragers (Fig. 4; see Fig. S4 for drone tissues). Expression of Am5- 
ht2oc was higher than Am5-ht2fi expression, e.g., up to ten times in 
the nervous system and ~200 times in the hypopharyngeal glands. 
Interestingly, in the brain and in the hypopharyngeal gland, the 
expression of the Am5-ht2uAIII transcript was lower than the full- 
length version, whereas the opposite held true for the Am5-ht2f5 
and Am5-ht2pAII transcripts (Fig. 4). All transcripts showed similar 
expression patterns in the CNS of drones compared with female 
honeybees (Fig. 4A+B; Fig. S4). 

Functional analyses of Am5-HT 2 receptors in HEK 293 
cells 

5-HT 2 receptors are known to induce an increase in [Ca 2+ ] ; . We 
generated HEK 293 cell lines constitutively expressing either 
Am5-HT 2o( or Am5-HT 2 p to examine whether these honeybee 
receptors also induced serotonin-dependent Ca 2+ signals. In a first 
series of experiments, we analyzed the ligand specificity of both 



receptors. Only serotonin (10 |jM) increased [Ca in Am5- 
HT 2a - and Am5-HT 2 p-expressing cells, whereas no signals were 
observed after the application of 10 |J,M tyramine, octopamine, or 
dopamine (Fig. 5). Non-transfected cells did not respond to any of 
these biogenic amines (Fig. 5). A third cell line expressing the 
alternatively spliced variant Am5-HT 2o( Aiii did not respond to the 
application of any biogenic amine, including serotonin. Therefore, 
this construct most likely does not code for a bona fide GPCR. 

For both full-length receptors, the concentration-response 
relationship for serotonin was examined with concentrations 
ranging from 1 nM to 300 uM. In cell lines expressing either 
Am5-HT 2o( or Am5-HT 2 p, the serotonin effect was concentration- 
dependent and saturable, resulting in sigmoidal concentration- 
response curves (Fig. 6A). The half-maximal effective concentra- 
tions (EC50) were similar at 25.7 nM and 32.5 nM for Am5-HT 2a 
and Am5-HT 2 p, respectively. Maximal increase in [Ca 2+ ]; was 
observed at serotonin concentrations of S 1 uM. 

To establish pharmacological profiles for both 5-HT 2 receptors, 
the effects of various 5-HT receptor agonists and antagonists were 
tested. Two agonists, namely 5-methoxytryptamine and 8- 
hydroxy-2-(di-n-propylamino)tetralin (8-OH-DPAT), caused an 
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Figure 2. Amino-acid sequence alignment of Am5-HT 2 „, Am5-HT 2 p, and orthologous receptors from Drosophila melanogaster (Dm5- 
HT 2 ; accession no. CAA57429) and Panulirus interruptus (Pan5-HT 2 |s ; AAS57919). Identical residues between all four receptors are shown as 
white letters against black, whereas conservatively substituted residues are shaded. Putative trans-membrane domains (TM1 -7) are indicated by gray 
bars. Putative consensus sites for post-translational modifications in Am5-HT 2a( and Am5-HT 2 p are indicated by open and filled symbols, respectively: 
N-glycosylation sites are shown as triangles, PKA phosphorylation sites as diamonds, PKC phosphorylation sites as circles, palmitoylation sites as 
squared boxes. Underlined letters represent the region within the CPL3 of Am5-HT 2o( used to raise specific antibodies (see Fig. S2). The amino acid 
position is indicated on the right. 
doi:1 0.1 371 /journal.pone.0082407.g002 



increase in [Ca 2+ ]; in both receptor-expressing cell lines (Fig. 6B). 
At the Am5-HT 2Q[ receptor, the EC50 values were 70 nM and 
55.9 u,M for 5-methoxytryptamine and 8-OH-DPAT, respective- 
ly. At the Am5-HT 2 p receptor, the EC 50 for 5-methoxytryptamine 
was 60.4 nM, and for 8-OH-DPAT, the EC 50 was 561.5 nM. 
Neither 5-methoxytryptamine nor 8-OH-DPAT had an effect on 
non-transfected HEK 293 cells. In contrast to the previously 
mentioned substances, 5carboxamidotryptamine and 2,5-di- 
methoxy-4-iodoamphetamine even affected [Ca "*]; in non-trans- 
fected cells and thus were not examined with the receptor- 
expressing cell lines. 



The potency of 5-HT-receptor antagonists to inhibit serotonin- 
induced Ca + responses in Am5-HT 2a - and Am5-HT 2 p-expressing 
cells was examined in another series of experiments. The 
calculated IC 50 and maximal inhibitory effect (I max ) of the 
substances are summarized in Table 1, and examples of 
concentration-response relationships for mianserin and clozapine 
are shown in Figures 6C+D. Most antagonists were able to inhibit 
the activation of the Am5-HT 2a receptor by 30-60%. At the Am5- 
HT 2a receptor, antagonists showed decreasing inhibitory potency 
with: SB200646 = mianserin > cyproheptadine > methiothepin 
> clozapine > methysergide. 
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Figure 3. Phylogenetic relationship between various 5-HT receptors from vertebrates, mollusks, and arthropods. Alignments were 
performed with BioEdit (version 7.0.5) by using the core amino-acid sequences lacking the variable regions of the N- and C-terminus and the third 
cytoplasmic loop. Maximum parsimony analysis was calculated with MEGA4. The receptor sequences followed by their accession numbers are listed 
in the order illustrated: Homo sapiens (Hs5-HT 1B , acc. no. NP_000854), Danio rerio (Dr5-HT 1B , AAI63698), Hs5-HT, D (NP_000855), Hs5-HT 1E 
(NP_000856), Hs5-HT 1F (NP_000857), Dr5-HT 1F (NP_001 076425), Hs5-HT 1A (NP_000515), Dr5-HT 1A (NP_001 1 39238), Aplysia californica (AcS-HT,, 
AAM46088, AAC28786), Panulirus interruptus (Pi5-HT„ AY528822), Apis mellifera (Am5-HT„ CBI75449), Drosophila melanogaster (Dm5-HT 1A , 
CAA77570), Dm5-HT 1B (CAA77571), Hs5-HT 5A (NP_076917), Dr5-HT 5A , (NP_001 1 19882), Am5-HT 7 (AM07671 7), Dm5-HT 7 (NP_524599), Aplysia kurodai 
(Ak5-HT 7 , ACQ90247), Hs5-HT 7A (NP_000863), Hs5-HT 7D (NP_062873), Hs5-HT 7B (NP_062874), Hs5-HT 4 (NP_001 035259), Ac5-HT 2 (Q16951), Ac5-HT, 
(Q16950), Hs5-HT 6 (NP_000862), Hs5-HT 2B (NP_000858), Dr5-HT 2B (ABI18978), Hs5-HT 2A (NP_000612), Dr5-HT 2A (CAQ15355), Hs5-HT 2C (NP_000859), 
Dr5-HT 2C (CAX14715), Am5-HT 2 p (FR727108), Dm5-HT 2p (NM_1 41 548/NMJ 69229), Pi5-HT 2p (AY550910), Am5-HT 2c( (FR727107), and Dm5-HT 2oi 
(CAA57429). The divergent D. melanogaster ninaE-encoded rhodopsin 1 (DmninaE, NM_079683) and D. melanogaster FMRFamide receptor (DmFR, 
AAF47700) were used as out-groups. The numbers at the nodes of the branches represent the percentage bootstrap support (2,000 replications) for 
each branch. 

doi:1 0.1 371 /journal.pone.0082407.g003 
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Figure 4. Tissue-specific expression patterns of Am5-ht2 receptor genes in pollen foragers determined by quantitative real-time 
PCR. Transcript levels were normalized to Amef-1x. (A) Relative levels of receptor-gene mRNAs in the brain (n = 1 1). (B) Relative levels of receptor- 
gene mRNAs in the ventral nerve cord (Am5-ht2a. and Am5-ht2zAIII: n = 1 1 ; Am5-ht2fi and Am5-ht2[lAII: n = 10) (C) Relative levels of receptor-gene 
mRNAs in the hypopharyngeal gland {Am5-ht2z: n = 6; Am5-ht2xAIII: n = 5; Am5-ht2fl: n = 7; Am5-ht2fiAII: n = 6). (D) Relative levels of receptor-gene 
mRNAs in the head salivary gland (Am5-ht2x and Am5-ht2uAIII: n = 5; Am5-ht2fl and Am5-ht2flAII: n = 7). (E) Relative levels of receptor-gene mRNAs in 
the thoracic salivary gland (each: n = 5). (F) Relative levels of receptor-gene mRNAs in Malpighian tubules {Am5-ht2x: n = 5; Am5-ht2zAIII: n = 6; Am5 
ht2/l: n = 8; Am5-ht2[lAII: n = 7). Groups that differed significantly from one another in relative mRNA levels are indicated with different letters above 
the box plots (p < 0.05, Bonferroni's multiple comparison test). 
doi:1 0.1 371 /journal.pone.0082407.g004 



Interstingly, SB-200646, methiothepin, and methysergide 
showed no effects at the Am5-HT 2 p receptor, whereas clozapine, 
cyproheptadine, ketanserin, and mianserin were strongly inhibi- 
tory with a maximal inhibitory effect of &77%. In comparison 
with Am5-HT 2ol , the rank order of inhibitory potency at Am5- 
HT 2 p was different with: cyproheptadine > ketanserin > 
mianserin > clozapine. 

Discussion 

The intense arborization of serotonergic neurons in the 
honeybee CNS suggests that serotonin controls the activity of 
many neuronal circuits [1 1]. In order to obtain further insight into 
serotonergic signaling in the honeybee, the molecular and 
functional identity of the 5-HT receptor repertoire of the 
honeybee has to be unraveled. Concomitantly with the release 
of the honeybee genomic sequence [37], several 5-HT receptor 
genes were annotated [10]. These in silico analyses, however, lack 
experimental data related to pharmacological properties, expres- 
sion patterns, and physiological functions of the identified 
receptors. So far, only a 5-HT 1A receptor [32] and a 5-HT 7 
receptor [33] have been thoroughly characterized. With the 



analysis of the two paralogous 5-HT 2 receptors, Am5-HT 2ot and 
Am5-HT 2 p, the current work provides an important step to 
complement our knowledge of the serotonergic receptor system in 
the honeybee. 

Genomic organization and gene prediction 

The coding sequences of Am5-ht2a and Am5-ht2fi are heavily 
spliced and consist of seven and six exons, respectively. Two splice 
sites are located in similar positions, one in ECL1 and a second in 
ECL2. These two splice sites are not only conserved in orthologous 
genes of the fruit fly (Fig. S3) but also represent the positions from 
where the splice variants originate (Fig. 1). 

A comparison of cloned and annotated Am5-ht2-receptor gene 
sequences reveals several discrepancies. Whereas the overall exon- 
intron structure is similar, the annotated Am5-ht2a sequence 
contains a shorter exon V (Fig. 1). Obviously, native honeybee 
spliceosomes not only use exactly conserved splice donor and 
acceptor sites, which, however, are central elements in gene 
prediction algorithms [43] . The situation for the Am5-ht2fi gene is 
even more complicated. Alignment of the cloned cDNA with the 
honeybee genome sequence showed that: (1) the Am5-ht2f5 gene 
contains extremely large introns (> 100,000 bp), which are hardly 
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Figure 5. Biogenic amine (10 jilvl each) modulation of intracellular Ca 2 concentrations in HEK 293 cells constitutively expressing 
Am5-HT 2 „ or Am5-HT 2 B receptors and in non-transfected HEK 293 cells. Cells were loaded with the Ca 2+ -sensitive dye Fluo-4 and changes in 
[Ca 2+ ]j were measured fluorimetrically. All values of HEK 293 cells stably expressing Am5-HT 2 p (black bars) and the corresponding non-transfected 
HEK 293 cells (gray bars) were normalized to the serotonin response of the Am5-HT 2 p-expressing cell line. The values obtained for the Am5-HT 2a cell 
line (white bar) and non-transfected HEK 293 cells (striped bar) were normalized analogously. Data represent the mean ± SD of octuplicate 
determinations from a representative experiment. The values for serotonin-stimulated Am5-HT 2a - and Am5-HT 2 p-expressing cells are significantly 
different from all other values (one-way ANOVA, followed by Dunnett's multiple comparison test, p < 0.05). 
doi:1 0.1 371 /journal.pone.0082407.g005 



recognized by gene prediction software [44] ; (2) the Am5-ht2ji gene 
is dispersed over three genomic contigs; (3) the ORF of Am5-ht2fS is 
oriented in either a 5' to 3' or 3' to 5' orientation compared with 
deposited contig sequences. This obviously prevents a computer- 
based assembly of the correct ORF. As exemplified by Am5-ht2f}, 
gene annotations provide important information, but molecular 
cloning of cDNA remains a prerequisite to determine the 
functional properties of a protein. 

Phylogenetic analysis of 5-HT 2 receptors 

Phylogenetic analyses of various 5-HT receptors revealed a 
separation of two main groups (Fig. 3). One group comprises 
receptors coupled to the cAMP pathway, i.e., 5-HT! (Gi-coupled) 
and 5-HT 7 receptors (G s -coupled). The other group comprises 5- 
HT 2 receptors preferentially coupled to G q proteins. Interestingly, 
deuterostomian and protostomian 5-HT receptors do not form 
distinct branches suggesting common ancestors for each receptor 
class. In arthropods, usually one or two receptor subtypes exist per 
5-HT receptor class. In D. melanogaster, the DmS-HT^ and Dm5- 
HT 1B receptors result from a recent duplication event [12]. So far, 
one 5-HT! receptor has been characterized in the honeybee 
(Am5-HT 1A , [32]). In addition, one 5-HT 7 receptor has been 
found in D. melanogaster and A. mellifera [14,33]. Here, we have 
added two functional 5-HT 2 receptors to the list of honeybee 5- 
HT receptors. A receptor orthologous to Am5-HT 2a has been 
characterized in D. melanogaster (CG1056; [13]). The D. melanogaster 
genome contains another gene (CG42796) potentially encoding an 
ortholog of Am5-HT 2 p. This hypothesis, however, remains to be 
established experimentally. 



Am5-ht2 genes give rise to several transcripts by 
alternative splicing 

Despite the apparent functional diversity afforded by the 
existence of multiple 5-HT 2 receptor subtypes, the existence of 
functional splice variants of these receptors adds to the potential 
repertoire of 5-HT receptors an additional mode of regulation. 
This phenomenon was observed for several 5-HT 2 receptor 
subtypes (Ascaris suum 5-HT 2 : [45]; 5-HT 2A : [46]; 5-HT 2C : 
[47,48]). Splice variants differ in distribution, ligand-binding 
properties and functional signaling (for a review, see: [49]). 

For both honeybee 5-HT 2 receptor genes, only one splice 
product contains information for a complete GPCR (Fig. 1). 
Shorter versions are generated by exon skipping. Interestingly, 
Am5-ht2c/.AIII and Am5-ht2pAII have lost the exon coding for a 
central part of the receptor including TM3 and TM4 (Fig. 1). 
Several 5-HT 2 receptor splice variants affecting the same region 
were also reported in mammalian genes [46,47]. A deletion 
covering the region encoding CPL2 and TM4 leads to a frameshift 
in the downstream sequence and a premature stop codon of the rat 
5-HT 2 c receptor [47]. Similarly truncated forms were also 
detected in cDNA samples from mouse and human brain [47]. 
A variant of the human 5-HT 2A receptor contains an 118 bp 
insertion which produces a frameshift and a premature stop codon 
[46]. Unfortunately, nothing is known about the function(s) of 
truncated 5-HT 2 receptors in mammals. 

When heterologously expressed, the Am5-HT 2a Ain receptor is 
synthesized as a glycosylated and non-glycosylated protein 
suggesting that it is processed in the Golgi apparatus (Fig. S2A). 
Attempts to monitor functional receptor activity from this 
construct were unsuccessful. This might be attributable to the 
absence, in TM3, of critical residues that interact with the amino 
group of biogenic amines [50]. Although the protein is not 
functional on its own, it might assemble with a full-length receptor 
and thereby modulate the binding and/or signal transduction 
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Figure 6. Concentration-dependent effects of agonists and antagonists on [Ca ]j in Am5-HT 2 „- and Am5-HT 2 | ( -expressing HEK 293 
cells. [Ca 2 *], is depicted as relative fluorescence given in percent. Data represent the mean ± SEM of eight replicates from one experiment 
representative of at least two similar experiments. (A) Am5-HT 23t - or Am5-HT 2 p-expressing cells were incubated with increasing concentrations of 
serotonin. Values were normalized to the response measured at the highest serotonin concentration (= 100%) for each cell line. (B) Cells expressing 
Am5-HT 2a or Am5-HT 2 p were incubated with increasing concentrations of the agonists 5-methoxytryptamine (5-MeOT) and 8-hydroxy-DPAT (8-OH- 
DPAT). Values were normalized to the response measured at the highest 5-MeOT concentration (= 100%) for each cell line. (C, D) Cells expressing 
Am5-HT 2a or Am5-HT 2 p were incubated with 50 nM serotonin and increasing concentrations of the antagonist mianserin (C) or clozapine (D). Values 
were normalized to the serotonin response in the absence of antagonist ( = 100%). 
doi:1 0.1 371 /journal.pone.0082407.g006 



Table 1. Antagonist profile of Apis mellifera 5-HT 2 receptors. 





Am5-HT 2a 


Am5-HT 2o 


Am5-HT 2 |) 


Am5-HT 2 p 


antagonist 


IC 50 [uM] 


maximal inhibition 


IC 50 I(iM] 


maximal inhibition 


clozapine 


1.7 


56% 


6.1 


95% 


cyproheptadine 


0.7 


64% 


0.2 


77% 


ketanserin 




no effect 


1.5 


89% 


methiothepin 


0.6 


61% 




no effect 


methysergide 


2.2 


45% 




no effect 


mianserin 


0.2 


51% 


2.4 


76% 


SB-200646 


0.1 


32% 




no effect 



IC 50 values (potency) and relative efficacy were calculated from concentration-response curves for each drug. Efficacy is given as the maximal inhibition (%) of 
fluorescence compared with 50 nM serotonin in the absence of antagonist and calculated by curve fitting. Values are means of the results of two independent 
experiments, each performed in octuplicates. 
doi:1 0.1 371 /journal.pone.0082407.t001 
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properties in a heterodimeric complex. Such modes of regulation 
have been observed for a Caenorhabditis elegans dopamine receptor 
variant lacking TM6 and TM7 [51] and for mutants of the human 
dopamine D2 receptor [52]. 

Deletion of exon II in Am5-ht2fiAII leads to a frameshift in the 
ORF, which is terminated by an early stop codon. Thus, ORF1 of 
Am5-ht2/} A II codes for a protein that shares, with full-length Am5- 
HT 2 p, the N-terminal part down to ECL1 but that thereafter 
contains a completely divergent sequence. A second ORF of this 
splice variant codes for a protein identical to the full-length 
receptor from TM5 down to the C-terminus (Fig. 1). Whether any 
of these proteins can be heterologously expressed and eventually 
modify the properties of the full-length receptor remains to be 
examined in a forthcoming study. 

Functional characterization of Am5-HT 2 receptors 

Am5-HT 2Q[ and Am5-HT 2 p were functionally expressed in 
HEK 293 cells. Like their mammalian orthologs, both receptors 
couple to the G q signaling pathway. These results support the 
concept that the signaling pathway is conserved for a given 
receptor class, whereas pharmacological profiles might differ. 

With EC 50 values of 25.7 nM and 32.5 nM, respectively, Am5- 
HT 2a and Am5-HT 2 p displayed similar potencies for serotonin as 
those also reported for 5-HT 2o( receptors from D. melanogaster 
(Dm5-HT 2o( , 16 nM; [13]) and Calliphora vicina (Cv5-HT 2a , 
24 nM; [36]). In addition to serotonin, 5-methoxytryptamine also 
activates honeybee 5-HT 2 receptors but with 2-fold higher EC 50 
values (Fig. 6B). In the blowfly C. vicina, nanomolar concentrations 
of 5-methoxytryptamine activate only Cv5-HT 2o( (EG 50 = 67 nM) 
but not Cv5-HT 7 [36]. However, at higher concentrations 
(EC50 = 55.9 uM and 560 nM, respectively), both honeybee 
receptors are also activated by 8-OH-DPAT, as is Cv5-HT 2a 
(EC50 = 62 uM). In conclusion, both 5-methoxytryptamine and 8- 
OH-DPAT can be considered as non-selective 5-HT receptor 
agonists in the honeybee, because 5methoxytryptamine also 
activates the Am5-HT 1A receptor [32], and because 8-OH-DPAT 
is a poor agonist at the Am5-HT 7 receptor [33]. 

In general, the pharmacological profiles of both Am5-HT 2 
receptors are different for antagonists, especially for SB-200646, 
methiothepin, and methysergide, which display high potency at 
Am5-HT 2a but do not block Am5-HT 2 p. To be considered as an 
Am5-HT 2a -specific antagonist, however, these substances need to 
be re-examined at the Am5-HTi and Am5-HT 7 receptors. 
Clozapine displays similar potencies for both 5-HT 2 receptors. 
In contrast, mianserin has a 24fold higher potency at the Am5- 
HT 2a receptor (Table 1) and thus is a promising candidate for an 
Am5-HT 2o( -specific antagonist. Interestingly, ketanserin seems to 
be a specific antagonist for Am5-HT 2 p. However, these substances 
also remain to be tested with AmS-HTj and Am5-HT 7 . 

We have demonstrated that many drugs are active at several 5- 
HT receptor subtypes of the honeybee ([32,33]; this study). Thus, 
when applying an agonist to elicit a physiological or behavioral 
"serotonin effect", the possible activation of multiple 5-HT 
receptors or even other aminergic receptors simultaneously should 
not be overlooked. 

Localization of Am5-HT 2 receptors 

The expression patterns of only two 5-HT receptors have been 
reported in the honeybee. The Am5-HT 1A receptor is restricted to 
neuronal tissues [32]. Transcripts of Am5-HT 7 have been detected 
in neural and peripheral tissues [33]. Based on our qPCR data, 
both Am5-ht2 receptor genes share the widespread expression 
pattern with Am5-ht7. 



Transcript levels for Am5-ht2oc were clearly higher than those for 
Am5-ht2fi, especially in the hypopharyngeal gland (Fig. 4C). The 
hypopharyngeal gland plays an important role in the control of 
eusocial behavior and displays 'functional flexibility' in secreting 
either royal jelly compounds or digestive enzymes depending on 
the age and task of the worker bee [53,54]. Queen caste 
development strongly depends on royal jelly [55,56]. A tempting 
speculation is that glandular activity is modulated by serotonin, 
and thus, signaling via Am5-HT 2o( is crucial for generating the 
female reproductive caste in honeybees. 

The Am5-ht2[l gene is expressed in the nervous system of worker 
bees and drones. In D. melanogaster, the orthologous receptor 
(CG42796) is preferentially expressed in the male CNS [57]. Thus, 
5-HT 2 p receptors might be involved in different functional circuits 
in these insects. 

With the present work, we extend the basic knowledge of the 
serotonergic system of the honeybee by adding two receptors, 
Am5-HT 2o( und Am5-HT 2 p, to the growing list of well-character- 
ized honeybee GPCRs. The detailed knowledge of the pharma- 
cological properties of a receptor is essential for designing and 
conducting targeted pharmacological experiments (e.g., [32]). In 
combination with expression data, knowledge of the pharmaco- 
logical profiles of these receptors should now facilitate such studies 
for individual 5-HT receptor subtypes in the honeybee. 

Supporting Information 

Figure SI Sequences of primers and TaqMan probes 
(including 5'- and 3'-modifications; see material and 
methods) used for qPCR assays and the expected length 
of the resulting amplicons. The positions of primers on the 
cDNAs are schematically shown by open arrows, TaqMan probes 
are indicated as asterisks. 
(PDF) 

Figure S2 Western blot and immunocytochemical anal- 
yses of Am5-HT 2a - and Am5-HT2 a Am- ex P r essing cell 
lines. Anti-Am5-HT 2a antibodies were raised against a fusion 
protein containing part of the third cytoplasmic loop (CPL3; 
amino acid Arg3 40 to Glu45 2 ; see Fig. 1). The cDNA fragment was 
amplified by PCR with specific primers (Table SI). The fragment 
was cloned into pMAL-c2X vector (New England Biolabs, 
Frankfurt, Germany). The fusion protein was over-expressed in 
E. coli BL2 1 CP and purified by amylose affinity-chromatography 
(New England Biolabs). In collaboration with the Nachwuchs- 
gruppe Antikorper-Technologien (University of Potsdam, Ger- 
many), the fusion protein was used to immunize mice and to raise 
monoclonal antibodies. A second fusion protein containing a His- 
tag attached to the same receptor fragment was expressed from 
pET-30a vector (Novagen, Darmstadt, Germany) and used for 
testing the specificity of the monoclonal antibodies. Membrane 
proteins (10 ug protein per lane) of human embryonic kidney cells 
(HEK 293) expressing Am5-HT 2a -HA and Am5-HT 2o , Am -HA 
receptors (see below) were isolated as previously described 
(Thamm et al., 2010). Proteins were separated by SDS 
polyacrylamide gel electrophoresis on 10% or 12% gels and 
transferred to polyvinylidene fluoride membranes (Roth, Karls- 
ruhe, Germany). These membranes were blocked with 5% (w/v) 
dry milk in Tris-buffered saline containing Tween 20 (TBS-T, 
10 raM Tris-HCl, pH 7.5, 150 mM NaCl, 0.01% Tween 20) for 
30 min at room temperature, incubated either with specific anti- 
HA antibodies (Anti-HA High Affinity, Roche, Penzberg, 
Germany; dilution 1:5,000) or with receptor-specific antibodies 
(dilution 1:100) in TBS-T, washed with TBS-T, and finally 
incubated with secondary antibodies (1:5,000, anti-rat-HRP; 
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American Qualex, La Mirada, USA; 1:200, anti-mouse Alexa568; 
Invitrogen) for 1 h. Signals were visualized by enhanced 
chemiluminescence. (A) Western blot analyses of membrane 
proteins (10 ug per lane) of non-transfected HEK 293 cells (nt) 
and HEK 293 cells expressing either full-length Am5-HT 2a -HA 
(full) or Am5-HT 2a Ain-HA (AIII) proteins. Both anti-HA (dilution 
1:5,000, left) and anti-Am5-HT 2Q[ (culture supernatant 1:10, right) 
antibodies recognize bands of identical size in protein preparations 
from transfected cells. No bands were detected in protein 
preparations from non-transfected cells. (B) Similar staining 
patterns were observed in receptor-expressing cell lines with both 
antibodies (anti-HA and anti-Am5-HT 2a ). Scale bar 40 |tm). 
Unfortunately, the anti-Am5-HT 2a antibody did not work with 
native tissue, neither on Western Blots nor on fixed tissue sections. 
This is probably due to the very low endogenous expression level 
of Am5-HT 2a , especially in nervous tissue. 
(PDF) 

Figure S3 Comparison of splice sites in 5-ht2 genes 
(coding regions) of Apis mellifera and Drosophila 
melanogaster. Splice donor and acceptor sites were identified 
using the Splign alignment tool (http:/ /www.ncbi.nlm.nih.gov/ 
sutils/splign/splign.cgi). In the mRNA sequence alignments 
(A+B), which were created using the program BioEdit version 
7.1.3.0, splice sites are indicated by +++ and highlighted in yellow 
and green for 5-ht2 genes of the honeybee and D. melanogaster, 
respectively. (A) Am5-ht2cn. is the orthologue of the Dm5-ht2a gene 
(CGI 056) with which it has three introns in common. (B) Am5- 
ht2f3 is the orthologue of the D. melanogaster CG42796 gene with 
which it has two introns in common. (C) A pictogram of a GPCR 
with its seven transmembrane segments (green bars) is displayed. 
Arrows point to the relative positions of splice sites in the primary 
structures of the receptors. For each exon, the last amino acid 
residue is indicated and numbered according to its position in the 
deduced amino acid sequence. One splice site is conserved in all 
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